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Ethanol  pressure-swirl  sprays  issuing  in  a  hot-diluted  oxidizer  coflow  with  different  temperature  and 
composition  were  studied.  The  bulk  coflow  temperature  was  varied  together  with  the  oxygen  volume 
fraction.  The  bulk  coflow  temperature  was  changed  from  1480  K  to  1225  K  and  the  oxygen  volume 
fraction  from  7.1%  to  10.1%.  The  liquid  mass  flow  rates  were  chosen  to  yield  spray  flames  with  nearly 
identical  Weber  number.  Laser  Doppler  anemometry,  phase  Doppler  anemometry  and  coherent  anti- 
Stokes  Raman  scattering  were  applied  in  the  spray  region  and  the  coflow  inlet.  The  current  measure¬ 
ments  provide  a  thorough  description  of  the  spray  structure,  droplet  dispersion  and  gas  temperature 
fields  as  well  as  a  comprehensive  database  useful  for  validation  of  numerical  models.  Visual  observations 
of  the  flame  luminescence  reveal  that  the  lift-off  height  depends  on  the  liquid  mass  flow  rates  as  well  as 
the  coflow  conditions.  The  lift-off  height  is  shown  to  increase  for  lower  coflow  temperatures  and  higher 
liquid  mass  flow  rates.  It  is  found  that  lift-off  behavior  depends  on  the  droplet  convective,  vaporization 
and  chemical  time  scales  prior  to  ignition.  Phase  Doppler  anemometry  results  indicate  that  the  droplet 
mean  size  and  velocity  distribution  close  to  the  atomizer  are  not  influenced  by  the  coflow  conditions. 
A  flame-front  develops  at  the  outer  region  of  the  spray  where  a  low  density  of  large  droplets  are  present. 
A  significant  number  of  peak  temperatures  samples  above  2000  K  is  observed  at  this  location.  Decrease  of 
the  coflow  temperature  leads  to  a  reduction  of  the  local  peak  temperatures.  Closer  to  the  center  axis,  the 
local  mixture  composition  becomes  increasingly  rich  and  the  heat-release  is  lower  than  in  the  outer 
region. 

©  2014  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Moderate  or  Intense  Low-Oxygen  Dilution  (MILD)  combustion 
has  prompted  great  interest  as  a  new  clean  combustion  concept. 
The  operation  principle  is  based  on  the  dilution  of  fuel  with  pre¬ 
heated  air  and  high  temperature  combustion  gas  recirculation 
yielding  an  overall  small  temperature  increase  due  to  the  heat 
capacity  augmentation  and,  consequently,  low  NOx  emissions 
(Wiinning  and  Wiinning,  1997;  Cavalieri  and  de  Joannon,  2004). 
This  unique  combustion  mode  has  completely  different  character¬ 
istics  from  a  conventional  operating  method  concerning  flamma¬ 
bility  limits,  combustion  stability,  combustion  noise,  and  ignition 
process  and  has  been  studied  both  in  industrial  furnaces  and  labo¬ 
ratory  scale  setups  for  various  gaseous  fuels.  Issues  of  interest  for 
the  studies  included  flame  chemistry  (Medwell  and  Dally,  2012; 
Medwell  et  al.,  2007;  Ihme  and  See,  2011),  flow  dynamics 
(Oldenhof  et  al.,  2011,  2013;  Dally  et  al.,  2002),  stabilization 
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mechanisms  (Oldenhof  et  al.,  2010),  emissions  formation  and 
quantification  of  the  thermal  field  uniformity  under  various  oper¬ 
ational  conditions  (Coelho  and  Peters,  2001;  Verissimo  et  al„ 
2011,  2013). 

A  review  of  earlier  research  show  that  MILD  combustion  for  liq¬ 
uids  fuels  has  received  little  attention  (Tsuji  et  al.,  2002;  Weber 
et  al.,  2005;  Wu  et  al„  2007;  Derudi  and  Rota,  2011).  Studies  have 
been  performed  mainly  in  industrial  and  semi-industrial  scale  fur¬ 
naces  and  focused  on  the  sustainability  of  MILD  regime  for  heavy 
fuel  oils  as  well  as  the  factors  resulting  in  furnace  performance 
deterioration.  Tsuji  et  al.  (2002)  investigated  the  applicability  of 
heavy  fuel  oils  in  industrial  furnaces  operating  in  MILD  conditions. 
The  oxygen  content  and  the  pre-heated  air  temperature  were  var¬ 
ied  systematically  (Y02  =  3-15%  and  T  =  573-1553  K)  to  investigate 
the  changing  state  of  the  flame  form  and  the  resulting  pollutant 
emission.  A  strong  dependence  of  the  flame  stability  and  chemilu¬ 
minescence  with  the  preheated  air  and  oxygen  dilution  was 
observed  as  well  as  a  substantial  increase  of  NOx  production  for 
heavy  fuel  oils.  Other  studies  (Weber  et  al.,  2005;  Wu  et  al., 
2007)  in  industrial  scale  test  furnaces  operating  in  MILD  conditions 
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showed  a  homogenized  temperature  field  in  furnaces.  However, 
only  for  certain  nozzle  configurations  and  fuel  oil  grades  a  high- 
temperature  process  without  the  penalty  of  increased  NOx  emis¬ 
sions  was  achievable.  Derudi  and  Rota  (201 1 )  carried  out  an  exper¬ 
imental  study  on  a  dual-nozzle  laboratory  scale  burner  arranged  in 
a  crossflow  configuration.  Pure  liquid  n-octane  and  other  mixtures 
involving  hydrocarbons  commonly  found  in  practical  fuels  were 
used.  The  experiments  were  designed  to  map  the  dependence  of 
the  MILD  combustion  region  with  two  main  operating  parameters: 
average  combustion  chamber  temperature  and  dilution  ratio. 
Results  on  NOx  and  CO  emissions  showed  that  MILD  conditions 
of  liquid  hydrocarbons  was  attainable  for  lower  combustion  gas 
recirculation  ratio  as  compared  to  the  LPG  fuel  counterparts.  Addi¬ 
tionally,  increase  of  the  oxygen  dilution  was  witnessed  to  have  a 
slight  positive  effect  on  pollutant  emissions  for  lower  average  com¬ 
bustion  chamber  temperature. 

The  aforementioned  studies  clearly  show  the  practical  impor¬ 
tance  of  the  technology,  however,  many  aspects  in  the  near-injec- 
tor  region  are  not  investigated.  In  the  spray  region,  the  fuel  stream 
first  has  to  be  shattered  into  an  ensemble  of  droplets  with  a  desir¬ 
able  distribution  of  sizes  and  velocities  in  order  to  achieve  the 
required  rate  of  vaporization,  chemical  heat-release,  levels  of  con¬ 
version  and  pollutant  emissions.  Instabilities  at  the  liquid-gas 
interface  in  the  near-atomizer  region  lead  to  primary  and  second¬ 
ary  break-up.  The  liquid  fragments  in  surrounding  gas  form  a  dis¬ 
persed  multiphase  flow,  where  droplets  are  evaporating  and 
subject  to  turbulent  dispersion.  The  dispersed  droplets  modify 
the  gas-phase  turbulence  and  the  released  fuel  vapor  depending 
on  the  relative  speed  of  evaporation,  mixing  and  combustion, 
burns  in  diffusion  flames  around  individual  droplets  or  clusters 
of  droplets,  or  in  partially  premixed  gaseous  flames  (Jenny  et  al., 
2012).  The  temperature  and  composition  of  the  oxidizer  influence 
the  time  scales  of  phase  change  and  reaction,  making  a  spray  flame 
in  conventional  combustion  different  from  that  of  a  hot-diluted 
combustion  regime. 

The  present  work  concerns  a  fundamental  experimental  study. 
The  objective  is  to  reach  understanding  of  the  spray  region  by 
using  a  laboratory  scale  burner  developed  to  study  spray  combus¬ 
tion  in  MILD  conditions.  The  burner  design  retains  the  relevant 
physical  processes  of  practical  combustion  systems.  However,  the 
composition  of  the  gases  entrained  into  the  spray  are  controlled 
by  a  secondary  burner  rather  than  the  furnace  aerodynamics.  The 
spray  flame  is  unconfined  and  specifically  designed  to  allow  for 
easy  access  of  state-of-the-art  laser  diagnostics  (Barlow,  2007). 

Goix  et  al.  (1994),  Cessou  and  Stepowski  (1996),  and  Stepowski 
et  al.  (1994)  investigated  stabilization  of  flame  structures  above  an 
airblast  injector  fed  with  liquid  methanol.  The  qualitative  structure 
of  the  reacting  spray  and  its  reaction  zone  were  investigated  by 
Mie  scattering  imaging  and  OH*  Planar  LIF.  Phase  Doppler  ane- 
mometry  was  used  to  measure  droplet  size  statistics.  The  results 
show  two  outer  reaction  zones  where  strong  OH  signals  are  pres¬ 
ent.  It  is  observed  that  flame  stabilization  occurs  in  the  large  scale 
mixing  structures  at  the  edge  of  the  jet  where  the  small  droplets 
have  sufficient  time  to  vaporize,  mix  with  entrained  air,  and  burn. 
A  discussion  was  presented  concerning  the  lift-off  position,  moti¬ 
vated  by  a  mixture-fraction  formulation.  In  an  effort  to  establish 
general  features  of  combusting  sprays  in  hot-diluted  coflow  with¬ 
out  the  influence  of  the  near-injector  dense  region,  O’Loughlin  and 
Masri  (2011)  used  a  nebuliser  instead  of  a  pressure  atomizer. 
Simultaneous  high-speed  OH-CH*  Planar  LIF  and  droplet  Mie  scat¬ 
tering  show  that  ignition  OH  kernels  formation  and  growth  is  the 
mechanism  of  spray  flame  stabilization  in  hot-diluted  conditions. 
The  oxygen  mole  fraction  in  the  coflow  (12%  by  volume)  of 
O’Loughlin  and  Masri  (2011)  is,  however,  rather  high  compared 
to  other  studies  of  combustion  systems  operating  in  MILD 
conditions  (Derudi  and  Rota,  2011;  Tsuji  et  al.,  2002). 


The  present  study  is  an  extension  of  earlier  work  by  the 
authors  on  the  structure  of  pressure-swirl  spray  flames  in  conven¬ 
tional  and  hot-diluted  combustion  regime  (Correia  Rodrigues 
et  al.,  2014).  The  burner  facility,  atomizer  type  and  the  liquid  fuel 
used  in  this  study  were  identical.  A  commercial  pressure-swirl 
atomizer  was  used  for  generating  the  spray.  In  this  way  the  liquid 
atomization  and  gas  flow  in  the  combustion  chamber  are  very 
strongly  coupled  and  experiments  remain  close  to  the  most 
common  configuration  used  in  industrial  practice.  All  measure¬ 
ments  are  made  up  to  locations  as  close  to  the  atomizer  as 
possible  to  unveil  features  of  the  near-injector  region.  Ethanol 
was  used  as  fuel  because  of  the  well-defined  physical  properties 
and  well-documented  chemical  mechanism  (Marinov,  1999; 
Rohl  and  Peters,  2009). 

In  the  authors’  previous  study  (Correia  Rodrigues  et  al.,  2014), 
measurements  of  profiles  of  droplet  size  and  velocity  and  gas- 
phase  temperature  were  performed  in  an  axisymmetric  spray  issu¬ 
ing  in  air  and  in  hot-diluted  coflow.  The  results  show  that  in  a 
reacting  spray  in  air  coflow,  an  inner  and  an  outer  flame-front 
are  present.  In  hot-diluted  coflow,  weakening  of  the  inner  flame- 
front  occurs  due  to  the  local  gaseous  mixture  richness.  At  the  loca¬ 
tion  of  the  outer  flame-front  temperature  samples  above  2000  I< 
persist.  In  this  study,  three  spray  issuing  in  hot-diluted  coflows 
with  different  temperature  and  composition  were  studied.  The 
present  study  complements  our  previous  work  (Correia 
Rodrigues  et  al.,  2014)  and  contributes  to  the  current  literature 
by:  (a)  describing  changes  in  the  flame  structure  and  identifying 
the  driving  parameter  for  peak  temperature  reduction  in  the  outer 
flame-front  region,  (b)  providing  a  comprehensive  set  of  measure¬ 
ments  that  extend  the  initial  database  for  model  validation  of  poly- 
dispersed  sprays,  and  (c)  provide  insight  into  the  physical 
mechanisms  controlling  the  spray  flame  lift-off  height  in  hot- 
diluted  conditions. 

This  paper  is  structured  as  follows:  Section  2  briefly  describes 
the  burner  facility  and  measurement  techniques  used  in  the  pres¬ 
ent  study.  Section  3  presents  the  input  parameters  as  well  as  the 
coflow  boundary  conditions  for  three  test  cases.  Section  4  presents 
the  results  on  spray  flames  visual  appearance,  lift-off  behavior, 
droplet  distribution  and  spray  flame  structure.  A  summary  of  the 
spray  phenomenology  and  a  discussion  of  the  spray  flame  stabil¬ 
ization  mechanisms  is  also  given.  In  Section  5,  conclusions  are 
presented. 


2.  Experimental  method 

2.1.  Burner  facility 

Experiments  are  carried  out  in  a  burner  facility  (Fig.  1)  at  Delft 
University  of  Technology.  It  consists  of  an  ethanol  spray  (C2H5OH) 
issuing  in  a  coflow  of  hot  combustion  products.  The  hot  combus¬ 
tion  products  are  generated  by  a  secondary  burner.  It  consists  of 
a  hexagonal  packing  arrangement  of  identical  5  mm  diameter  ver¬ 
tical  pipes  with  each  end  welded  to  perforated  plates.  At  half¬ 
height  of  each  individual  vertical  pipe,  four  0.5  mm  holes  were 
drilled  in  a  cross-configuration  (see  section  A-A  in  Fig.  1).  The 
incoming  air  flows  through  the  vertical  pipes  and  mixes  with 
Dutch  Natural  Gas  (DNG)  forced  through  the  0.5  mm  holes.  The 
mixture  then  forms  a  matrix  of  small  flames  stabilized  on  the  per¬ 
forated  plate.  A  vertical  pipe  with  270  mm  length  and  160  mm 
inner  diameter  prevents  the  coflow  of  hot  combustion  products 
from  mixing  with  the  surrounding  air.  As  the  hot  combustion  prod¬ 
ucts  flow  through  the  vertical  pipe  their  temperature  is  decreasing 
due  to  the  heat  loss  to  the  surroundings  by  radiation  from  the 
vertical  pipe.  Hereafter,  the  coflow  of  hot  combustion  products  will 
be  refer  as  ‘hot-diluted  coflow’. 
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Pressure-swirl  atomizer 


TOP  VIEW 


Fig.  1.  Schematic  of  the  burner  facility  with  relevant  dimensions  in  millimeters. 


The  ethanol  spray  is  produced  by  a  commercial  pressure-swirl 
solid  cone  atomizer  (Delavan  WDA  0.5  GPH)  with  an  0.21  mm  exit 
orifice  and  a  60°  spray  angle  mounted  on  a  2  mm  inner  diameter 
stainless  steel  pipe  through  which  liquid  fuel  is  fed  to  the  atomizer. 
During  operation  the  burner  components  absorb  a  substantial 
amount  of  heat  due  to  the  high  coflow  temperature.  In  order  to 
avoid  pre-vaporization  of  the  ethanol,  a  water-cooling  system  con¬ 
sisting  of  a  helical  copper  coil  was  built  around  the  stainless  steel 
pipe  through  which  ethanol  fuel  is  fed.  Ethanol  temperature  was 
monitored  near  the  atomizer  nozzle  thread  by  a  thermocouple 
Thermocoax  Type  K.  The  liquid  injection  pressure  was  measured 
online  with  a  pressure  transducer  Omegadyne  MMA.  Acquisition 
of  pressure  data  required  dedicated  electronics  hardware  to  con¬ 
nect  the  pressure  transducer  output  to  the  TSI  E1C  external  data 
collection  system  board. 

The  air  and  DNG  mass  flow  rates  are  set  by  Bronkhorst  mass 
flow  controllers  with  ±0.5%  inaccuracy  of  full  scale.  The  mass  flow 
rate  of  the  liquid  fuel  was  set  by  a  Coriolis  flow  meter  with  ±0.5% 
inaccuracy.  The  burner  was  mounted  on  a  three-axis  orthogonal 
traverse  system  controlled  by  stepper  motors.  Experimental  equip¬ 
ment  used  remained  in  place  and  are  described  in  Section  2.2. 
Throughout  this  paper,  a  r,  z-coordinate  system  is  used  with  its 
origin  at  the  atomizer  orifice  exit  plane.  The  convention  of  aligning 
‘z’  with  the  axial  direction  and  Y  with  the  transverse  direction  is 
adopted. 

2.2.  Measurement  techniques 

In  this  study,  five  measurements  techniques  were  separately 
employed:  laser  Doppler  anemometry,  phase  Doppler  anemome- 
try,  coherent  anti-Stokes  Raman  spectroscopy  and  flue  gas  ana¬ 
lyser.  A  brief  description  of  each  system  used  is  given  below. 

2.2.1.  Laser  Doppler  Anemometry  (LDA) 

LDA  measurements  were  performed  at  z  =  0  mm,  i.e.  at  the 
atomizer  orifice  exit  plane  and,  also,  along  the  spray  axis  where 


no  droplets  are  present.  The  results  provide  the  input  boundary 
conditions  needed  for  numerical  simulations  as  well  as  relevant 
information  to  understand  of  the  spray  flame  behavior.  A  two- 
dimensional  LDA  system  was  employed  to  measure  simulta¬ 
neously  the  axial  and  radial  velocity  components  Uz  and  Ur, 
respectively.  The  LDA  system  consists  of  a  1 0  W  Argon-ion  laser, 
Bragg  cell,  probe  head  and  a  processor  unit.  A  TSI  350  mm  focal 
length  lens  collimated  the  laser  beams  emerging  from  the  probe 
head  yielding  a  probe  volume  of  0.17  mm  in  diameter  and 
2.40  mm  in  length.  The  scattered  light  from  the  seeding  particles 
in  the  probe  volume  was  collected  in  backscattering  mode  and 
guided  to  a  FSA-3000  process  unit.  An  additional  TSI  EIC  external 
data  collection  system  was  used  to  record  the  analog  output  of 
the  pressure  transducer. 

The  air  flow  was  seeded  with  alumina  particles  (A1203)  by 
means  of  a  cyclone  system.  The  nominal  particle  median  diameter 
is  of  1  pm.  At  each  measurement  point  104  samples  were  acquired 
to  ensure  adequate  convergence  of  the  statistics.  Velocity  statistics 
were  computed  using  the  transit  time  as  a  weighting  factor  to 
compensate  for  the  velocity  bias  (Albrecht  et  al.,  2002).  Conse¬ 
quently,  the  velocity  statistics  should  be  interpreted  as  Favre-aver- 
aged  quantities  (Oldenhof  et  al.,  2011).  According  to  the  results  of 
the  detailed  analysis  of  the  various  error  sources  in  LDA  measure¬ 
ments  as  reported  by  Absil  (1995)  we  assumed  that  the  mean 
velocity  components  are  accurate  up  to  2%  of  the  maximum  flow 
velocity.  The  Reynolds  normal  stresses  and  the  Reynolds  shear 
stress  are  accurate  up  to  5%  and  8%  of  their  local  maximum  values, 
respectively. 

2.2.2.  Phase  Doppler  Anemometry  (PDA) 

A  two-dimensional  PDA  system  was  employed  to  measure 
simultaneous  droplet  velocity  and  size  statistics  in  the  spray 
region.  The  laser,  Bragg  cell  and  probe  head  used  were  the  same 
as  in  the  LDA  system  described  above.  The  receiver  was  set  up  in 
first-order  refraction  and  had  a  650  mm  focal  length  receiving  lens. 
The  measurement  volume  was  truncated  by  a  150  pm  slit  and  the 


312 


H.  Correia  Rodrigues  et  al.  / International  Journal  of  Heat  and  Fluid  Flow  51  (2015)  309-323 


signal  fed  into  a  FSA-4000  processor  unit.  The  processor  unit  was 
also  connected  to  a  TSI  EIC  external  data  collection  to  collect  the 
analog  output  of  the  pressure  transducer. 

Velocity  and  size  statistics  of  a  polydispersed  ensemble  of  drop¬ 
lets  can  be  significantly  affected  by  adjustable  parameters  of  the 
PDA  measurement  system  (Albrecht  et  al.,  2002).  To  ensure  the 
reproducibility  of  the  measurements,  a  parametric  study  of  the 
adjustable  operational  parameters  that  influence  the  droplet  diam¬ 
eter,  droplet  number  density  and  velocity  were  performed  prior  to 
the  experimental  campaigns.  The  selected  settings  were  kept  con¬ 
stant  throughout  the  measurements  campaigns.  The  droplet  den¬ 
sity  varies  throughout  the  spray  domain  and  the  number  of 
samples  collected  should  be  sufficient  to  achieve  convergence  of 
the  statistical  moments  within  a  feasible  measurement  time.  A 
total  number  of  104  samples  were  collected  at  positions  with  data 
rates  lower  than  102  Hz.  For  all  the  remain  ensembles  of  3  x  104 
samples  were  acquired. 

Based  on  the  work  of  Widmann  and  Presser  (2002),  the  inaccu¬ 
racy  associated  with  the  droplet  mean  velocity  is  about  4%.  Oscil¬ 
lations  in  the  PDA  response  are  found  to  occur  for  small  droplets 
and  results  in  different  sizing  error  depending  on  the  droplet  diam¬ 
eter.  According  to  the  work  of  Taylor  et  al.  (1994),  sizing  errors  for 
droplets  below  1  pm  are  on  the  order  of  60%.  Above  1  pm,  the  siz¬ 
ing  errors  are  approximately  8%.  For  further  details  on  PDA  errors 
and  uncertainties  we  refer  to  the  work  of  Friedmann  and 
Renksizbulut  (1999).  The  droplet  diameter  is  estimated  based  on 
the  phase  shift  between  the  Doppler  signals  and  the  droplet  refrac¬ 
tive  index.  For  all  cases  studied,  an  estimation  of  the  droplet  tem¬ 
perature  based  on  theoretical  models  of  single-component  droplet 
vaporization  (Birouk  and  Gbkalp,  2006)  indicates  that  droplets  are 
at  boiling  temperature  at  the  first  measurement  station.  Thus,  the 
effect  of  temperature  on  the  refractive  index  of  ethanol  is  negligi¬ 
ble  and  no  systematic  error  in  the  droplet  diameter  is  expected. 

Smaller  droplets  follow  faithfully  the  gas  motion  and  were  used 
to  determine  the  gas-phase  statistics.  The  selection  of  the  droplet 
size  representative  of  the  gas  motion  was  estimated  by  computing 
the  first-and  second-order  statistics  with  different  droplet  size 
thresholds.  For  all  cases  in  this  study,  droplets  with  diameter 
below  6  pm  yielded  minimal  differences  in  the  first  and  second 
order  statistics  and  where  chosen  as  representative  of  the  gas 
motion.  As  will  be  shown  later  in  Section  4.8,  estimates  of  the 
Stokes  number  for  different  droplet  size  classes  additionally  con¬ 
firm  the  selected  droplet  size  threshold  as  suitable  for  tracers  of 
the  gaseous  phase. 

2.2.3.  Coherent  Anti-Stokes  Raman  Scattering  (CARS) 

The  CARS  system  was  employed  to  evaluate  the  gas-phase  tem¬ 
perature  statistics  at  z  =  0  mm  and  in  the  spray  region.  A  brief 
description  of  the  CARS  system  and  data  post-processing  proce¬ 
dure  is  reported  here  but  for  further  details  the  reader  is  referred 
to  van  Veen  and  Roekaerts  (2005,  2003).  The  ro-vibrational  CARS 
system  is  based  on  an  injection-seeded,  frequency-doubled 
Nd:YAG  laser  (Spectron  SL805  SLM),  which  yields  500  mj  per  pulse 
at  532  nm  with  a  pulse  duration  of  12  ns  at  10  Hz  repetition  rate.  In 
a  planar-boxcars  phase-matching  configuration,  a  probe  volume  of 
700  pm  length  and  35  pm  diameter  is  obtained.  The  generated 
CARS  radiation  at  the  probe  volume  is  focused  onto  the  entrance 
of  an  echelle  spectrometer,  and  the  spectra  dispersed  on  a  CCD 
detector  with  1100  x  330  pixels.  The  spectra  are  contained  in 
two  strips  of  1100  intensity  values,  which  are  digitized  by  an 
18-bit  AD  converter  and  stored.  Dacapo  software  is  used  to  nor¬ 
malize  the  CARS  spectra  using  a  nonresonant  spectrum  and  fit  to 
a  library  of  theoretical,  temperature-dependent  spectra.  Previous 
experiments  show  that  single-shot  imprecision  of  the  system  is 
1-4%  over  a  range  from  2000  K  to  300  K.  The  inaccuracy  is 
estimated  to  be  20  I<  (van  Veen  and  Roekaerts,  2005,  2003). 


Three  considerations  were  taken  into  account  in  the  CARS  mea¬ 
surements.  First,  the  incidental  presence  of  a  droplet  along  the 
beams  paths  manifests  itself  as  an  increase  of  the  N2  CARS  resonant 
signal  background.  A  code  was  implemented  in  Dacapo  software  to 
evaluate  and  substract  a  linear  spectral  background  profile,  inter¬ 
polated  between  values  found  at  the  left-  and  right-hand  sides  of 
the  spectrum  before  proceeding  to  the  temperature  evaluation. 
Secondly,  a  steep  temperature  gradient  may  occur  over  the  probe 
volume,  particularly  near  the  reaction  zone  of  the  flame.  In  the 
presence  of  such  a  gradient,  a  substantially  lower  temperature 
than  the  mean  temperature  will  be  calculated  due  to  the  density 
effect.  This  was  corrected  based  on  a  weighted  fit  according  to 
the  procedure  suggested  by  van  Veen  and  Roekaerts  (2003).  Last, 
at  certain  locations  of  the  spray  region,  nonresonant  signals  were 
observed  in  addition  to  resonant  N2  signal.  Based  on  the  species 
present,  only  ethanol  -  having  a  large  nonresonant  susceptibility 
-  can  account  for  this  signal.  A  quantitative  analysis  is  not  possible, 
but  by  the  mere  presence  of  nonresonant  signals  regions  with  fuel 
vapor  can  be  identified  qualitatively. 

2.2.4.  Flue  gas  analyser 

Oxygen  volume  fraction  measurements  were  performed  at 
z  =  0  mm  with  a  Testo  335  fluegas  analyser  having  a  specified  inac¬ 
curacy  of  ±0.20%. 

3.  Cases  description 

The  operating  conditions  as  well  as  relevant  dimensionless 
numbers  of  each  individual  flame  studied  are  given  in  Table  1. 
The  positions  were  the  measurement  techniques  were  employed 
are  summarized  in  Table  2.  Three  spray  flames  in  different  coflow 
conditions  were  thoroughly  studied.  The  three  flames  are  denoted 
H!,  Hn  and  H„[.  In  addition,  some  results  are  reported  for  two  aux¬ 
iliary  cases  Hfux  and  H*ux  corresponding  to  sprays  with  higher 
liquid  mass  flow  rates  in  identical  hot-coflow  conditions  as  ^ 
and  Hu,  respectively.  These  are  included  to  elucidate  the  droplet 
velocity  statistics  in  the  near  atomizer  region  and  help  the  discus¬ 
sion  of  the  flame  lift-off  behavior  presented  in  Section  4.9.  Only 
PDA  measurements  were  made  for  these  two  cases.  UCf,  Fcf  and 
(Xo2)  represent  the  mass  average  velocity,  temperature  and  oxygen 
volume  fraction  of  the  coflow,  respectively.  The  values  were  deter¬ 
mined  by  the  measured  velocity,  temperature  and  oxygen  volume 
fraction  profiles  at  z  =  0  mm,  i.e.  at  the  atomizer  orifice  exit  plane. 

The  Weber  number  is  defined  as  pgUjetl  j a  and  is  identical  for  the 
three  spray  flames  studied.  I  is  the  characteristic  liquid  film  thick¬ 
ness  and  it  was  assumed  to  be  equal  to  the  nozzle  diameter,  i.e. 
210  pm.  Ujet  represents  the  difference  between  the  estimated 
liquid  jet  velocity  from  the  liquid  mass  flow  rate  and  Ucf. 

The  coflow  velocity  profiles  at  z  =  0  mm  are  displayed  on  Fig.  2. 
On  the  left-hand  side,  the  coflow  axial  mean  velocity  (Uz)  and  the 
respective  r.m.s.  fluctuations  u'z  are  presented.  The  same  quantities 
are  shown  for  the  radial  velocity  component  on  the  right-hand 
side.  It  is  observed  that,  on  opposite  sides  of  the  spray,  for  all  cases 
the  magnitude  of  ( Ur )  sharply  increases  until  a  maxima  in  the  near¬ 
wall  of  the  central  rod  and  then  gradually  decreases  outwards.  This 
indicates  that  the  hot-diluted  coflow  has  a  strong  inward  velocity 
component  and  is  being  convected  towards  the  spray  axis  filling 
in  the  void  space  between  the  ligaments  and  droplets  formed  at 
the  center. 

The  r.m.s.  fluctuations  profiles  of  the  axial  and  radial  velocity 
components,  u'z  and  u'r,  show  a  sharp  increase  between  r=  70  and 
80  mm.  This  indicates  that  a  shear  layer  is  formed  between  the 
hot-diluted  coflow  and  the  laboratory  air.  Thus,  the  gradual  mixing 
of  laboratory  air  with  the  hot-diluted  coflow  causes  changes  in  the 
radial  profiles  of  temperature  and  oxygen  concentration  and, 
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Table  1 

Burner  operating  conditions  and  dimensionless  numbers  of  the  flames  studied.  The  pressure  fluctuations  are  less  than  0.2%  for  all  flames. 


Case 

Coflow 

Spray 

&  i-i 

mair  [kg/h] 

riiDNG  [kg/h] 

He/ [m/s] 

Fc/IK] 

(x0!)  m 

mLiq  [kg/h] 

P  [bar] 

TUq  [K] 

Ujet  [m/s] 

We  [-] 

H, 

70 

3.28 

3.42 

1480 

7.1 

1.36 

12.0 

305 

14.0 

0.3 

0.035 

j_jAux 

1.53 

15.9 

304 

15.7 

0.4 

0.037 

Hu 

51 

2.13 

2.23 

1300 

9.3 

1.46 

11.5 

301 

14.9 

0.4 

0.048 

ijAux 

HII 

1.63 

14.6 

301 

16.7 

0.7 

0.051 

H,„ 

45 

1.75 

1.81 

1225 

10.1 

1.48 

11.7 

298 

15.1 

0.5 

0.052 

Table  2 

Summary  of  measurements  locations  for  cases  H,,  Hn  and  H„|. 


Measurement 

technique 

z  (mm) 

Comment 

LDA 

0,  20,  25,  30,  35,  40,  45, 
50 

Centerline,  coflow  radial 
profiles 

CARS 

0,  15,  20,  30,  40,  50,  60 

Coflow  radial  profiles,  spray 
region 

Flue  gas  analyser 

0 

Coflow  radial  profiles 

PDA 

8,  10,  12,  15,  20,  30,  35, 
40,  45 

Spray  region,  steps  of  1  mm 

ultimately,  perturbs  the  spray  flame.  It  will  be  shown  in  Section  4.6 
that  up  to  z  =  60  mm,  the  burner  pipe  diameter  is  wide  enough  to 
isolate  the  spray  flame  from  the  laboratory  air  entrainment. 

Fig.  3  displays  the  coflow  mean  temperature  (T)  and  r.m.s.  fluc¬ 
tuations  T'  as  well  as  the  oxygen  volume  fraction  (X02)  for  the  three 
cases  presented  in  Table  1.  For  the  cases  studied,  (X02)  radial  pro¬ 
files  are  fairly  uniform  over  a  significant  radial  distance  in  the 
coflow,  in  the  temperature  profiles,  asymmetries  are  seen  that 
exceed  the  inaccuracy  of  the  CARS  system.  The  maximum  differ¬ 
ence  between  the  { T)max  at  opposite  sides  of  the  spray  is  observed 
in  case  Hu  and  amounts  to  34  K. 


4.  Results  and  discussion 

4.3.  Spray  flames  visual  appearance  and  behavior 

To  show  the  overall  flame  appearance,  color  pictures  are  dis¬ 
played  in  Fig.  4.  The  flames  were  photographed  using  a  digital  cam¬ 
era  (Canon  EOS  X5)  and  have  0.6  s  exposure  time  and  identical 
aperture.  The  right  inset  shows  a  close-up  of  the  spray  region 


-80  -60  -40  -20  0  20  40  60  80 

Radial  position  [mm] 


and  initial  development  of  the  flame.  For  all  the  flames  presented, 
the  field-of-view  of  the  inset  measures  170  x  160  mm2 
(width  x  height). 

It  is  observed  that  the  atomizer  yields  conical  spray  flames. 
Based  on  the  visual  appearance,  the  spray  flames  in  hot-diluted 
coflow  are  comprised  of  three  regions:  jet  liquid  break-up  with 
no  visible  flame  luminescence,  a  bluish  reaction  zone  and  a  yellow 
zone  farther  downstream.  The  images  shows  that  flame  Hm  has  a 
stronger  luminescence.  Furthermore,  a  shorter  flame  length  is 
observed  as  compared  to  flames  Hi  and  Hn.  For  all  the  cases,  the 
insets  show  a  region  of  bright  luminescence  in  the  outer  regions 
of  the  spray. 

A  study  was  performed  to  determine  the  lift-off  height  depen¬ 
dence  on  the  liquid  mass  flow  rates  for  the  coflow  conditions  of 
case  Hn  and  Hm.  Experiments  were  also  conducted  for  spray 
flames  in  air  coflow  and  have  been  reported  in  Correia  Rodrigues 
et  al.  (2014).  In  the  current  study,  the  procedure  to  determine 
the  lift-off  height  was  based  on  average  flame  luminescence  from 
the  photographs.  It  is  defined  as  the  axial  distance  from  the  atom¬ 
izer  orifice  exit  plane,  z  =  0  mm,  to  the  lowest  axial  position  of  the 
bluish  reaction  zone  along  the  spray  axis. 

The  outcome  is  presented  in  Fig.  5  and  the  coflow  conditions  are 
presented  in  the  inset.  The  test  cases  of  this  study  are  identified  in 
the  figure  by  the  dashed  squares.  Case  An,  a  reacting  spray  in  air 
coflow,  was  previously  described  by  the  authors  (Correia 
Rodrigues  et  al.,  2014)  and  is  also  presented.  When  the  liquid  mass 
flow  rates  were  below  1.10  kg/h,  pre-vaporization  of  the  ethanol 
occured  yielding  a  yellowish  flame  anchored  at  the  atomizer  orifice 
exit  plane.  Liquid  mass  flow  rates  higher  than  1.70  kg/h  could  not 
be  achieved  with  the  current  experimental  setup. 

As  shown  in  Fig.  5,  the  lift-off  height  of  sprays  in  hot-diluted 
conditions  depends  linearly  on  the  liquid  mass  flow  rate  whereas 
in  air  coflow  a  non-linear  dependence  is  observed.  This  is  evidence 


Radial  position  [mm] 


Fig.  2.  Coflow  mean  velocities  and  r.m.s.  fluctuations  profiles  of  the  axial  and  radial  velocity  components  at  z  =  0  mm.  ( 


:  Hm). 
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Fig.  3.  Mean  oxygen  volume  fraction  :.X0l),  mean  temperature  (T)  and  root-mean- 
square  temperature  fluctuations  T  profiles  at  z  =  0  mm.  ( — □ — :  Hi,  — O — •  Hu, 
HIU). 

for  the  existence  of  different  stabilization  mechanics  of  the  spray 
flames  in  hot-diluted  coflow.  Furthermore,  it  is  observed  that  in 
the  coflow  conditions  of  case  Hm,  the  lift-off  height  is  more  sensi¬ 
tive  to  the  liquid  mass  flow  rate  than  in  the  conditions  of  cases  Hi 
and  Hn. 

Secondary  breakup  is  not  expected  to  play  an  important  role  in 
these  flames.  For  the  sprays  in  air  coflow  the  Weber  number  varies 
from  2  to  3.4  and  does  not  exceed  a  critical  Weber  number  of  12 
(based  on  Wierzba  (1990)).  Thus,  in  air  coflow  the  secondary 
breakup  is  expected  to  be  absent.  In  hot-diluted  conditions,  Weber 
numbers  are  smaller  and  within  the  range  of  liquid  mass  flow  rates 
studied  only  vary  from  0.2  to  0.5.  Although,  to  our  knowledge  in 
the  literature  no  information  is  available  on  the  secondary  breakup 
modes  in  the  presence  of  strong  evaporation  rates,  the  ratio  of  the 
aerodynamic  forces  and  surface  tension  remains  significantly 
below  unity  and,  therefore,  no  changes  on  the  secondary  breakup 
regime  are  expected. 


Liquid  mass  flow  rate  [kg/hr] 


Fig.  5.  Lift-off  height  dependence  on  the  liquid  fuel  mass  flow  rate  based  on  visual 
observations.  The  bulk  properties  of  the  coflow  inlet  are  indicated  in  the  table. 

4.2.  Droplet  size  and  mean  velocity  field 

Fig.  6  shows  droplet  Sauter  mean  diameter  (d32)  and  mean  axial 
(Uz)  and  radial  (Ur)  velocity  profiles  for  cases  Hb  Flu  and  Hm  at  sev¬ 
eral  axial  stations.  For  flame  Hi,  no  droplets  are  observed  above  the 
axial  station  of  30  mm.  Regions  in  the  spray  core  where  no  data  are 
plotted  correspond  to  positions  were  acquisition  of  sufficient 
amount  of  data  within  a  reasonable  measuring  time  was  not 
possible. 

The  experimental  results  indicate  that  the  initial  d32  and  veloc¬ 
ity  initial  profiles  as  well  as  subsequent  evolution  downstream  are 
similar  for  the  three  flames  studied.  At  z=  10  mm  the  sprays  are 
non-uniform  with  smaller  d32  near  the  centerline  and  large  d32  at 
the  periphery.  It  is  observed  that  droplets  have  a  zero  (Ur)  close 
to  the  centerline.  From  the  centerline  until  the  spray  periphery, 
(Ur)  magnitude  gradually  increases  until  a  maximum.  This  indi¬ 
cates  that  droplets  are  spreading  into  the  hot-diluted  coflow  and 
it  is  in  agreement  with  the  visual  appearance  of  the  spray  flames 
(see  Fig.  4  insets).  With  increasing  radial  distance,  (Uz)  increases 
until  a  maximum  and  then  decreases  to  nearly  the  same  magni¬ 
tude  as  in  the  center  region.  For  all  the  cases,  a  good  symmetry 
of  both  velocity  components  and  d32  is  observed  on  opposite  sides 
of  the  sprays  axis. 

With  increasing  axial  distance  from  the  atomizer,  the  center 
region  of  the  sprays  become  void  of  droplets.  As  expected,  this 


(a)  Flame  Hi 


(b)  Flame  Hu 


(c)  Flame  Hm 


Ucf  —  3.42  m/s 
Tcf  =  1480K 
<Xo2>  =  7.1% 


Ucf—  2.24  m/s 
Tct  =  1300  K 
(Xq2)  =  9.3% 


Ucf—  1-81  m/s 
Tcf  =  1225  K 
(Xq2)  =  10.1% 


Fig.  4.  Images  of  spray  flames  studied  with  0.6  s  exposure  time  and  identical  aperture.  The  bulk  properties  of  the  coflow  inlet  for  each  spray  flame  are  listed  below. 
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Fig.  6.  Droplets  Sauter  mean  diameter  d32  and  mean  axial  (Uz)  and  radial  (Ur)  velocity  profiles  at  several  axial  stations.  The  different  colored  vertical  axis  scales  display  the 
aforementioned  physical  quantities.  The  upper  and  lower  limits  have  been  kept  the  same  for  the  different  axial  stations.  (-  -  -  (14),  (Ur),  Q  d32). 


effect  is  more  pronounced  for  flame  Hi  since  Tc f  is  the  highest  and 
droplets  evaporate  faster.  A  gradual  reduction  of  the  droplet  veloc¬ 
ity  components  is  observed  due  to  the  momentum  exchange  of  the 
droplets  with  the  gaseous  phase.  An  increase  in  the  d32  between 
z  =  1 0  mm  and  z  =  30  mm  is  observed.  The  reason  for  is  that  small 
drag/inertial  effects  of  large  droplets  leads  to  the  displacement  of 
the  droplets  outwards,  progressively  increasing  the  liquid  fuel 
mass  available  in  the  spray  peripherical  regions  of  the  spray. 


4.3.  Spray  flames  initial  droplet  distribution 

To  demonstrate  the  importance  of  the  initial  conditions  on  the 
overall  spray  morphology  and  development,  the  results  of  polydis- 
perse  ensemble  of  droplets  at  z=  10  mm  is  plotted  in  Fig.  7.  The 
three  radial  positions  selected  correspond  to  the  spray  centerline, 
the  radial  position  where  the  droplet  number  concentration  is 
maximum  (see  Fig.  8)  and  the  spray  edge.  For  case  Hu  no  droplets 
are  present  at  r  =  0  mm  and  the  radial  position  of  2  mm  was 
selected  as  representative  of  the  centerline.  Droplet  size  frequency 
distribution  curves  (lefty-axis)  are  overlaped  with  the  correspond¬ 
ing  cumulative  droplet  distribution  curves  (right  y-axis).  The  distri¬ 
bution  curves  based  on  number  of  droplets  per  size  class  are 
represented  by  black  lines  whereas  gray  lines  refer  to  volume  size 
and  can  be  related  to  mass  by  the  liquid  fuel  density.  Note  that  the 
droplet  distribution  based  on  volume  are  skewed  to  the  right  as 
expected  due  to  the  weighting  effect  of  the  larger  droplets. 

The  results  show  clear  differences  between  the  droplet  distribu¬ 
tions  at  the  radial  positions  presented,  however,  the  trend  among 
the  three  cases  studied  are  in  common.  Flame  H!  shows  that,  at 
r  =  2  mm,  50%  of  the  droplets  are  below  5  pm  indicating  that  a  high 
droplet  surface  area  is  available.  50%  of  the  droplet  mass  is  con¬ 
tained  in  the  fraction  below  18  pm.  At  larger  radial  distances  from 
the  spray  centerline,  r  =  6mm  and  r=13mm,  the  spread  of  the 
droplet  number  distribution  increases  and  more  larger  droplets 


are  observed.  The  droplet  diameter  below  which  50%  of  the  mass 
is  contained  also  increases. 

Regarding  flames  Hu  and  Hm,  it  is  observed  that  at  all  the  three 
radial  positions  larger  droplets  are  more  abundant  in  the  droplet 
size  distributions  than  in  H^  At  the  centerline,  r  =  0  mm,  50%  of 
the  droplet  mass  resides  in  droplets  below  28  pm  and  24  pm  for 
flame  Hn  and  Hnl,  respectively.  At  r  =  6  mm  it  can  be  noticed  that, 
as  compared  to  flame  Hi,  droplet  number  cumulative  distributions 
do  not  change  significantly.  For  the  radial  position  of  r=  13  mm 
and  14  mm  a  considerable  number  of  large  droplets  is  clearly 
observed  and,  therefore,  more  fuel  mass  is  available  at  the  spray 
periphery  for  cases  Hn  and  Hm. 

4.4.  Droplet  number  density 

The  collection  of  liquid  droplets  spreads  through  the  gaseous 
environment  according  to  their  initial  momentum  and  droplet- 
gas  interactions.  Radial  profiles  of  the  droplet  number  density 
(nd)  at  several  axial  stations  are  shown  in  Fig.  8  in  order  to  eluci¬ 
date  the  degree  of  dispersion  and  identify  regions  where  gas-phase 
composition  and  momentum  may  be  significantly  influenced  by 
the  droplets  presence. 

For  all  cases,  at  z=  10  mm,  nd  is  nonuniform  across  the  radial 
traverse  and  droplets  mainly  are  concentrated  at  r  =  6  mm.  The 
subsequent  development  shows  a  broadening  of  the  spray  with 
increasing  axial  distance  and  a  significant  reduction  of  the  droplet 
number  density.  Furthermore,  in  the  region  close  to  the  spray  axis 
no  droplets  are  present.  This  effect  is  more  significant  for  case  Hi.  In 
all  the  cases  studied,  the  largest  radial  distance  from  the  spray  axis 
at  which  droplets  are  observed  is  nearly  the  same. 

4.5.  Gas-phase  flow  field 

The  evolution  of  the  gas-phase  mean  axial  (Uz)g  and  radial  (Ur)g 
velocity  components  and  normal  ( uzuz)g  and  shear  ( uzur)g  Reynolds 
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Fig.  7.  Frequency  distribution  curves  (solid  lines)  and  cumulative  droplet  size  distribution  (dashed  lines)  based  on  number  (d)  and  volume  (d3)  at  z  =  10  mm  for  flames  HIt  H„ 
and  Hm. 


stress  profiles  (Hanjalic  and  Launder,  201 1 )  between  four  axial  sta¬ 
tions  is  displayed  in  Fig.  9.  The  subscript  “g”  denotes  for  the  gas- 
phase.  Each  point  presented  in  Fig.  9  is  based  on  an  ensemble  of, 
at  least,  100  samples.  This  threshold  was  found  by  investigating 
the  convergence  of  the  velocity  mean  and  higher  order  statistical 
moments  of  droplets  smaller  than  6  pm  as  the  number  of  samples 
increases.  Notice  that  profiles  of  the  physical  quantities  presented 
are  similar  for  the  different  cases  and  the  curves  have  been  shifted 
for  clarity.  Due  to  the  striking  similarity  among  the  profiles  the  fol¬ 
lowing  observations  are  common  for  all  the  three  cases  studied. 

The  { Uz)g  profiles  indicate  that  at  z  =  10  mm  a  peak  value  occurs 
at  about  r=  6  mm.  In  the  region  around  the  centerline,  (Uz)g  mag¬ 
nitude  is  close  to  the  small  droplets  axial  velocity  (see  Fig.  6)  indi¬ 
cating  that  droplets  move  along  with  the  gas  flow  field  and  no 
substantial  momentum  exchange  through  drag  takes  place  in  the 
center.  At  downstream  axial  stations,  the  (Uz)  magnitude  gradually 
decreases  across  the  spray  domain.  This  occurs  as  a  result  of  the 
turbulent  mixing  with  the  entrained  hot-diluted  coflow.  Moreover, 
the  number  droplet  density  is  relatively  high  at  this  radial  position. 
Hence,  a  modulation  of  the  gas  velocity  field  by  the  droplets  is 
expected. 

The  ( Ur)g  magnitude  is  one  order  of  magnitude  smaller  than 
(Uz)g  at  z=  10  mm.  Further  downstream,  a  reduction  of  the  radial 
velocity  component  takes  place  to  a  nearly  nil  value.  This  shows 
that  the  gas  flow  field  move  mainly  upwards  and  droplets  shoot 


Fig.  8.  Variation  of  droplet  number  density  nd  with  the  radial  distance  measured  at 
different  axial  positions  for  flames  HIt  Hn  and  Hin. 
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Fig.  9.  Gas-phase  mean  axial  (Uz)g  and  radial  velocity  (Ur)g  along  with  Reynolds-stresses  (uzur)g  and  (uzuz)g  at  four  axial  stations  for  each  flame  studied.  (- 
— O — :  z  =  1 5  mm,  — 0 — :  z  =  20  mm, - 1 — :  z  =  30  mm). 


z=  10  mm. 


outwards.  Depending  on  their  initial  size  diameter  and  momen¬ 
tum,  the  droplet  trajectory  path  is  deflected  towards  the  direction 
parallel  to  the  spray  axis. 

Both  Reynolds-stresses,  { uzuz)g  and  { uzur)g ,  are  high  close  to  the 
exit  plane  of  the  atomizer.  Moreover,  a  large  difference  between 
{ uzuz)g  and  (uzur)g  components  is  observed,  with  (uzuz)g  substan¬ 
tially  higher  than  { uzur)g .  That  can  be  attributed  to  the  effect  of 
liquid  break-up  in  the  dense  region.  Additionally,  it  indicates  that 
significant  mixing  takes  place  in  the  near  atomizer  region.  A  fast 
decay  of  fluctuation  levels  until  30  mm  is  observed  resulting  in 
small  anisotropy  in  the  regions  farther  downstream. 

Since  the  number  of  droplets  is  very  low  in  the  center  region  of 
the  spray,  LDA  measurements  were  conducted  along  the  spray  cen¬ 
terline  to  determine  the  gas-phase  velocity.  The  LDA  system  can 
measure  the  velocities  of  both  ‘seeding’  particles  that  faithfully  fol¬ 
low  the  gas-phase  and  of  droplets.  Therefore,  ‘cross-talk’  from  the 
velocity  distribution  of  the  droplets  to  that  of  the  particles  poten¬ 
tially  results  in  velocity  bias  in  the  gas-phase  statistics  (Hardalupas 
and  Taylor,  2004).  The  degree  of  bias  will  depend  on  the  correlation 
between  the  particle  and  droplet  arrival  rate.  In  view  of  that,  mea¬ 
surements  were  performed  only  in  positions  were  no  droplets  are 
present,  i.e.  in  the  center  region  void  of  droplets.  Additionally,  prior 
to  any  further  estimation  of  the  flow  velocity  component,  velocity 
distribution  were  meticulously  studied  to  look  for  the  evidence  of 
outliers. 

The  evolution  of  gas-phase  axial  mean  { Uz)g  and  r.m.s.  ufl(Uz)g 
fluctuations  along  the  spray  axis,  is  illustrated  in  Fig.  10  for  the 
three  flames  studied.  For  all  the  flames,  { Uz)g  gradually  increases 
until  z  =  60  mm.  As  the  flow  is  converted  downstream,  a  mean 
gas  acceleration  in  the  center  regions  takes  place  due  to  turbulent 
diffusion  of  momentum  from  regions  away  from  the  center.  The 
peak  values  of  { Uz)g  for  cases  Hn  and  Hln  exhibit  a  somewhat  smal¬ 
ler  magnitude  than  case  FL,  reflecting  the  slower  mixing  rates  in 
the  center.  The  results  of  uf /(Uz)g  indicate  that  turbulence  levels 
in  the  centerline  are  substantially  higher  for  flame  H,. 


4.6.  Gas-phase  temperature  field 

The  radial  profiles  of  the  gas-phase  temperature  statistics  at 
four  axial  locations  are  shown  in  Fig.  11.  On  the  left-hand  side, 
the  gas  mean  temperature  (T)  is  presented  along  with  the  corre¬ 
sponding  mean  temperature  gradient  (d(T)/dr).  On  the  right-hand 
side,  the  r.m.s.  fluctuations  T  are  presented. 

At  z  =  15  mm,  the  (T)  profiles  show  a  rise  from  the  coflow  to  a 
single  peak  centered  around  r  =  14  mm  for  all  cases.  From  the  peak 
temperature  location  towards  the  spray  axis,  (T)  gradually 


Fig.  10.  Gas-phase  mean  axial  velocity  ( Uz)g  and  normalized  r.m.s  fluctuations 
uf/(llz)g  along  the  spray  axis  for  flames  Hi,  Hu  and  Hm.  ( — □ — :  H,,  — O — :  Hu, 
— O— :  Hm). 
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decreases  to  values  lower  than  the  ones  observed  in  the  coflow. 
However,  in  the  region  close  to  the  spray  centerline,  (T)  is  signifi¬ 
cantly  higher  than  the  ethanol  ignition  temperature  (rjg  =  653  K). 
The  observed  { T )  profile  shape  can  be  attributed  to  the  effects  of 
evaporation  and  mixing  prior  to  ignition  as  follows:  due  to  the 
liquid  atomization,  an  inward  convective  motion  of  the  coflow  is 
generated.  The  hot-diluted  coflow  is  entrained  into  the  spray  and 
provides  heat.  The  fuel  vapor  arises  at  the  cost  of  heat  losses  from 
the  oxidizing  environment  yielding  lower  gas-phase  temperatures 
towards  the  spray  axis.  The  mixture  composition  evolves  from  lean 
at  the  outer  parts  of  the  spray  towards  more  rich  at  the  center. 
Since  the  gas  has  relatively  high  temperature,  the  lean  mixture  at 
the  spray  outer  parts  ignites  and  the  heat-release  yields  a  peak  in 
the  (T).  As  one  approaches  the  spray  axis,  the  mean  heat-release 
is  lower  and  thus  a  decrease  of  (T)  takes  place.  The  peak  (T)  located 
at  r=  14  mm  is  nearly  identical  for  the  three  cases  studied.  This 
results  from  the  combined  effects  of  the  coflow  composition  and 
temperature  on  the  flame  adiabatic  temperature.  The  increase  of 
the  oxygen  dilution  in  the  oxidizer  stream  yields  a  reduction  of 
the  stoichiometric  adiabatic  temperature.  However,  a  higher 


oxygen  dilution  in  the  coflow  is  accompanied  by  an  increase  of 
the  overall  temperature  of  the  oxidizer  stream.  Thus,  the  two 
effects  cancel  out  and  no  substantial  differences  among  the  three 
cases  are  observed. 

At  z  =  20  mm,  three  observations  can  be  made.  First,  for  all 
cases,  the  peak  (T)  shifts  towards  larger  radial  distance  and  it  is 
located  in  a  region  of  large  droplets  and  low  droplet  number  den¬ 
sity  (see  Fig.  8,  17  <  r<  20  mm).  This  effect  results  in  a  decrease  of 
the  mean  temperature  gradient  from  the  peak  (T)  towards  the  cen¬ 
ter  as  compared  to  z  =  10  mm.  Secondly,  the  mean  temperatures  in 
the  region  0  <  r<  15  mm  are  highest  for  case  Hi.  This  observation  is 
in  agreement  with  the  sharper  reduction  of  the  droplet  concentra¬ 
tion  in  this  region  for  case  Hi  (see  Fig.  8).  Lastly,  throughout  the 
flame  region,  the  r.rn.s.  temperature  fluctuation  is  rather  similar 
among  the  three  cases  and  the  magnitude  is  overall  higher  by 
75  K  than  in  the  hot-diluted  coflow. 

At  higher  axial  stations,  z  =  30  and  40  mm,  the  displacement  of 
the  peak  (T)  towards  larger  radial  distances  proceeds  due  to  the 
continuous  spreading  of  the  large  droplets  into  the  hot-diluted 
coflow.  However,  the  distance  between  the  peak  and  the  spray  axis 
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Fig.  11.  Gas-phase  mean  temperatures  (T)  (left)  and  r.rn.s.  fluctuations  T  (right)  at  four  axial  stations.  The  dashed  gray  lines  correspond  to  the  corresponding  mean 
temperature  gradients. 
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is  smaller  for  case  Hi.  This  indicates  that  spray  flame  Hi  is  narrower 
than  the  sprays  flames  in  Hn  or  Hm.  This  is  in  agreement  with  the 
photographs  presented  in  Fig.  4. 

At  z  =  40  mm,  a  peculiarity  in  the  radial  profiles  of  (T)  and  T  can 
be  observed.  This  occurs  more  prominently  in  cases  Hn  and  Hm  and, 
therefore,  a  more  detailed  description  is  given.  For  these  cases, 
between  the  radial  positions  r  =  0  and  r=10mm,  V  increases  to 
a  local  peak  at  r  =  10  mm  and  then  decreases  to  smaller  values  until 
r=  15  mm.  Also,  a  significant  change  in  the  (T)  slope  takes  place  at 
r  =  10  mm.  An  authors’  early  study  describes  the  spray  flame  struc¬ 
ture  in  air  (case  An)  and  in  hot-diluted  coflow  (case  Hn).  In  air 
coflow,  an  inner  and  an  outer  flame-front  are  present  starting  from 
the  lift-off  position  and  diverging  with  increasing  downstream 
location  (Correia  Rodrigues  et  al.,  2014).  In  the  case  of  hot-diluted 
coflow  (case  Hn),  the  inner  flame-front  is  weak.  Due  to  the  local 
mixture  richness  the  heat-release  is  significantly  smaller  as  com¬ 
pared  to  the  outer  flame-front  and  a  local  maxima  in  (T)  is  not 
observed.  The  results  of  Fig.  1 1  show  that  for  different  coflow  tem¬ 
perature  and  oxygen  dilution  an  inner  flame-front  exists  in  all 
cases.  It  would  be  worthwhile  to  investigate  this  in  more  detail 
using  other  laser  diagnostics  such  as  OH-Planar  Laser  Induced 
Fluorescence. 

Fig.  1 1  shows  additional  information  concerning  the  influence 
of  air  entrainment  on  the  coflow  inlet  temperature  profiles.  It 
can  be  noticed  that  the  coflow  temperature  at  r  =  40  mm  remains 
nearly  the  same  throughout  the  axial  stations  presented.  This 
shows  that,  for  the  range  of  radial  distances  displayed,  the  burner 
pipe  diameter  is  wide  enough  to  isolate  the  spray  flame  from  the 
laboratory  air  entrainment.  CARS  measurements  at  higher  axial 
stations  (not  shown)  indicate  that  above  z  =  60mm  the  outer 
region  of  the  spray  flame  is  perturbated  by  the  laboratory  air.  Thus, 
below  z  =  60  mm,  the  spray  flame  is  stabilized  only  by  the  hot- 
diluted  coflow.  Above  the  axial  station  of  60  mm,  the  oxidizer 
stream  consists  of  a  mixture  of  laboratory  air  together  with  hot- 
diluted  coflow  and  reaction  zone  weakening  effects  may  be  present 
(Medwell  et  al.,  2009). 

4.7.  Gas-phase  temperature  histograms 

The  normalized  temperature  histograms  of  cases  Hi,  Hu  and  Hm 
at  z  =  40mm  are  presented  in  Fig.  12.  Five  radial  positions  were 
selected  for  each  case.  The  histogram  corresponding  to  the  radial 
position  in  the  coflow  core  region  is  represented  by  a  thicker  line. 
The  remaining  positions  correspond  to  the  spray  region.  A  bin  size 
of  50  K  was  selected  for  all  the  histograms  presented.  The  total 
number  of  samples  in  each  histogram  is  1000,  except  for 
r=  0  mm  of  case  H,  where  the  number  of  valid  samples  was  120. 
For  each  of  the  five  radial  positions  selected,  the  corresponding 
mean  (T),  root-mean-square  temperature  fluctuation  V  as  well  as 
skewness  (S)  and  kurtosis  (K)  of  the  temperature  distribution  are 
listed  in  the  table  inset. 

At  the  radial  positions  where  the  peak  (T)  is  observed  (r  =  20,  24 
and  24  mm  for  cases  Hi,  Hu  and  Hm,  respectively),  the  temperature 
distribution  are  nearly  symmetric  (S  =  0)  and  broad  with  tempera¬ 
tures  ranging  from  2400  K  down  to  1400  K.  A  considerable  amount 
of  samples  are  over  2000  K  for  all  cases.  It  is  observed  that  the 
number  of  valid  samples  above  2000  K  is  greater  for  case  Hi.  This 
is  related  to  the  higher  temperature  of  the  oxidizer  coflow. 

For  all  cases,  the  radial  position  of  15  mm  is  located  between 
the  inner  and  outer  flame-front.  As  expected,  lower  temperature 
samples  are  present  in  the  temperature  distributions.  As  compared 
to  the  radial  positions  where  the  peak  ( T )  occurs,  a  substantial 
reduction  of  the  number  of  valid  samples  above  2000  K  is 
observed.  For  case  H^  a  substantial  amount  of  valid  samples  above 
2000  K  still  persist  at  this  position,  yet,  it  is  significantly  reduced  as 
compared  to  r  =  20  mm.  These  results  suggest  that  the  coflow 
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Fig.  12.  Normalized  temperature  histograms  for  five  radial  positions  at  z  =  40  mm 
for  cases  Hi,  Hn  and  Hm. 

temperature  is  the  driving  parameter  with  respect  to  the  outer 
flame-front  temperature  peaks. 

At  the  spray  axis,  r  =  0  mm,  the  temperature  distributions  are 
symmetric  and  samples  are  below  1200  K  for  all  cases.  Tempera¬ 
ture  histograms  at  z  =  20  mm  and  30  mm  (not  shown)  have  virtu¬ 
ally  identical  r.m.s.  and  skewness  value.  This  indicates  that  in  the 
region  close  to  the  spray  axis  only  turbulent  mixing  occurs  and 
the  gas-phase  is  rather  homogeneous. 


4.8.  Droplet-size  class  properties 

Figs.  13  and  14  show  the  mean  axial  (Uz)  and  radial  (Ur)  velocity 
components  of  four  droplet  size-classes  at  three  axial  stations, 
z=  10,  20  and  30  mm.  The  class  width  was  chosen  as  10  pm.  All 
the  points  are  an  average  of  at  least  200  samples.  The  droplets  in 
the  size  class  between  0  and  6  pm  are  considered  to  be  flow 
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Fig.  13.  Radial  profiles  of  the  conditioned  droplet  mean  axial  velocity  ( Uz )  for  three  axial  stations.  ( 
20  pm  sj  d  <  30  pm,  — 0 — :  30  pm  <  d  <  40  pm, - 1 — :  40  pm  <  d  <  50  pm). 
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Fig.  14.  Radial  profiles  of  the  conditioned  droplet  mean  axial  velocity  (Uz)  for  three 
20  pm  <  d  <  30  pm,  — 0 — :  30  pm  <  d  <  40  pm, - 1 — :  40  pm  ^  d  <  50  pm). 

tracers.  The  mean  velocity  of  this  size  class  thus  represents  the  gas- 
phase  mean  velocity.  The  difference  between  the  mean  velocity  for 
a  particular  size  class  and  the  gas-phase  mean  velocity  will  be 
referred  to  as  slip  velocity. 

The  development  of  droplet  (Uz)  along  the  several  axial  stations 
shows  similar  trends  for  cases  Hp  Hn  and  Hm.  At  z=  10  mm,  an 
increase  of  droplet  slip  velocity  with  the  size-class  is  observed.  Far¬ 
ther  downstream,  slip  velocity  gradually  decreases  for  the  different 
droplet  size-classes.  Identical  observations  can  be  made  regarding 


axial  stations.  ( — 0 — :  Gas-phase  (d .  6  pm),  — □ — :  10  pm  '  d  <  20  pm,  — Q — : 

the  droplet  (Ur)  component.  This  occurs  because  the  momentum  of 
the  liquid  droplets  is  transferred  to  the  gaseous  phase  through  drag 
and  evaporation. 

4.9.  Discussion  of  the  flame  structure  and  lift-off  behavior 

This  section  presents  a  discussion  on  the  droplet  dispersion, 
flame  structure  and  lift-off  behavior  of  spray  flame  in  hot-diluted 
coflow  conditions.  Before  proceeding  with  the  analysis,  it  is  of 
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Table  3 

Relevant  dimensionless  numbers  and  relaxation  time  scales  for  different  droplet  size  classes. 


Size-class  (urn) 

Cases 

H, 

H„ 

H,„ 

iiAux 

Hl 

ljAux 

HII 

Wed  [-] 

10-20 

0.0 

0.0 

0.0 

0.0 

0.0 

20-30 

0.1 

0.1 

0.1 

0.1 

0.1 

30-40 

0.1 

0.2 

0.2 

0.2 

0.2 

40-50 

0.2 

0.2 

0.2 

0.3 

0.3 

Red  [-1 

10-20 

0 

1 

1 

1 

1 

20-30 

1 

2 

2 

2 

2 

30-40 

2 

3 

4 

3 

4 

40-50 

3 

4 

5 

4 

5 

Std  1-1 

0-6 

0.5 

0.5 

0.5 

0.5 

0.5 

10-20 

11 

13 

13 

14 

14 

20-30 

32 

36 

38 

36 

40 

30-40 

62 

70 

75 

70 

78 

40-50 

102 

116 

124 

115 

129 

id  [ms] 

0-6 

0.01 

0.01 

0.01 

0.01 

0.01 

10-20 

0.17 

0.18 

0.19 

0.17 

0.18 

20-30 

0.47 

0.50 

0.53 

0.47 

0.50 

30-40 

0.93 

0.99 

1.63 

0.93 

0.99 

40-50 

1.54 

1.63 

1.72 

1.54 

1.63 

interest  to  discuss  in  detail  relevant  dimensionless  numbers  char¬ 
acterizing  the  droplet  motion. 

4.9.1.  Relevant  dimensionless  numbers 

Some  relevant  dimensionless  numbers  are  listed  in  Table  3  for 
each  of  the  size  classes  and  for  the  droplets  smaller  than  6  pm  that 
are  representative  of  the  statistics  of  the  gas-phase.  The  droplet 
Weber  number  (Wed)  was  determined  as  pgU2dc/rr.  Here,  pg  repre¬ 
sents  the  gaseous  phase  density,  dc  is  the  mean  droplet  diameter  of 
the  selected  size  class  and  \Jj  the  difference  between  the  maximum 
droplet  (Uz)  measured  at  z=  10  mm  and  the  bulk  coflow  velocity. 
The  gas  properties  were  determined  from  the  coflow  bulk  temper¬ 
ature  and  composition.  As  can  be  observed  in  Table  3,  Wed  is  well 
below  unity  for  all  droplet  size  classes  suggesting  that  small  drop¬ 
lets  resulting  from  the  secondary  breakup  are  absent  for  all  five 
cases. 

The  droplet  Reynolds  number  (Red)  was  determined  as  Ujdc/vg. 
Assuming  that  droplet  are  spherical  and  move  in  a  viscous  flow  at  a 
very  low  Reynolds  number,  the  relaxation  time  (id)  can  be  deter¬ 
mined  as  pddd  / 18  pg.  Here,  pd  and  p„  are  the  droplet  density  and 
the  gaseous  phase  dynamic  viscosity,  respectively.  As  shown  in 
Table  3,  for  droplets  below  10  pm  the  Red  is  well  below  one 
whereas  for  droplets  bigger  than  30  pm  the  highest  value  encoun¬ 
tered  for  Red  is  of  5.  Hence,  larger  droplets  incur  an  underestima¬ 
tion  of  the  droplet  relaxation  time  by  at  most  7%.  On  the  other 
hand,  droplets  with  a  diameter  smaller  than  10  pm  satisfy  the 
assumption  of  Stokes  flow. 

The  Stokes  number  (Std)  is  defined  as  the  ratio  of  the  character¬ 
istic  time  of  a  droplet  to  a  characteristic  time  of  the  flow  and  indi¬ 
cates  the  different  regimes  of  droplet  dispersion  for  each  droplet 
size  class  (Eaton  and  Fessler,  1994;  Wood  et  al.,  2005).  The  flow 
characteristic  time  was  defined  as  da,/Uat  where  da,  is  the  atomizer 
nozzle  diameter  and  Uat  the  estimated  liquid  jet  velocity  in  the 
atomizer  nozzle.  Table  3  shows  that  for  all  cases  studied,  droplets 
bigger  than  10  pm  have  Std»  1  and  therefore  do  not  follow 
changes  in  the  flow  field.  For  droplets  with  diameter  smaller  than 
6  pm,  the  droplet  response  time  is  faster  than  the  flow  character¬ 
istic  time  scale  and,  therefore,  behave  as  flow  tracers. 

4.9.2.  Flame  structure 

Fig.  15  presents  a  schematic  of  the  spray  structure  in  hot- 
diluted  coflow  conditions  and,  together  with  Table  3,  can  be  used 
to  describe  the  droplet  dispersion  and  the  flame  structure.  As  a 


result  of  the  spray  breakup,  droplets  spread  outwards  and  an  outer 
flame-front  develops  wherein  low  droplet  concentration  and  large 
droplets  prevail.  The  outward  ballistic  motion  of  the  droplets 
results  from  the  droplet  initial  radial  velocity  as  well  as  the  corre¬ 
spondent  large  St  values  and  is  the  mechanism  responsible  for  the 
flame  inverted  conical  shape.  As  shown  in  Table  3,  droplets  bigger 
than  10  pm  have  Stokes  number  sufficiently  large  to  not  follow 
changes  in  the  flow  field.  Thus,  as  the  spray  evolves  downstream, 
large  droplets  move  outwards  and  a  flammable  mixture  resulting 
from  the  entrained  hot  oxidizing  coflow  and  the  fuel  vapor  is 
formed  at  increasing  larger  radial  distances.  A  comparison  of  cases 
Hi  and  Hm  additionally  shows  highest  Stokes  number  for  droplets 
within  the  same  size  class  when  the  coflow  temperature  is  lowest 
(Hm).  Hence,  lower  coflow  temperatures  yields  a  deeper  penetra¬ 
tion  of  droplets  into  the  hot-diluted  coflow  and,  consequently, 
the  flames  become  wider.  This  is  in  agreement  with  the  visual 
observations,  namely  the  flame  of  case  Hm  is  slightly  wider  than 
that  of  Hi  (see  Fig.  4  insets). 

4.9.3.  Characteristic  scales  and  lift-off  behavior 

The  trends  presented  in  Fig.  5  show  that  the  lift-off  height 
decreases  for  increasing  bulk  coflow  temperatures  while  it 
increases  with  increasing  liquid  mass  flow  rates.  This  section  pre¬ 
sents  a  possible  explanation  of  the  lift-off  behavior  of  spray  flames 
in  hot-diluted  coflow.  As  proposed  by  Delabroy  et  al.  (1998)  and 
Martinez-Ruiz  et  al.  (2013),  in  a  reacting  spray,  the  time  needed 
for  combustion  to  take  place  is  determined  by  the  time  scales  of 
three  phenomena:  (1)  droplet  vaporization,  (2)  droplet  convection, 
and  (3)  the  chemical  reaction.  The  comparison  of  the  different  time 
scales  allows  to  characterize  the  lift-off  height  behavior  of  the 
spray  flame.  Estimates  of  characteristic  numbers  for  different  drop¬ 
lets  size  classes  are  presented  in  Table  4.  The  following  discussion 
is  based  on  single  droplet  theory  and  the  interdroplet  interaction 
effects  are  neglected  (Zoby  et  al.,  2011;  Miyasaka  and  Law, 
1981).  For  all  the  calculations,  the  gas  properties  were  determined 
from  the  coflow  bulk  temperature  and  composition. 

The  droplet  vaporization  time  scale  ( tvap )  was  estimated  by  the 
d2  law  (Law,  1982;  Birouk  and  Gokalp,  2006).  It  is  considered  that 
droplets  emerge  from  the  atomizer  nozzle  at  the  boiling  point  and, 
therefore,  a  droplet  heat-up  period  is  absent.  For  an  evaporating 
spherical  droplet  in  a  flow,  tvap  can  be  written  as  dl/K.  K  is  the 
vaporization  constant  and  is  given  by  -4/ Pi[pmDmShm  ln(l  +  BM)]. 
This  model  accounts  for  the  slip  velocity  between  the  gas  and 
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Fig.  15.  Schematic  of  the  gas  and  droplet  motion  of  pressure-swirl  sprays  in  hot- 
diluted  coflow  conditions. 

droplet  motion.  The  gas-phase  properties  values  ‘m’  were  evalu¬ 
ated  by  the  1/3-rule  (Chen  and  Pereira,  1996).  D  is  the  binary  dif¬ 
fusion  coefficient  of  ethanol  and  air  at  condition  ( Tm ,  Ym ).  Sh  is 
the  Sherwood  number  and  BM  represents  the  Spalding  mass  trans¬ 
fer  number.  Immediately  after  the  atomizer  nozzle  exit  plane,  the 
droplets  travel  a  certain  path  and  leave  behind  a  wake  of  fuel  vapor 
that  depends  on  the  vaporization  rates.  Since  the  chemical  reaction 
cannot  begin  until  the  gaseous  fuel  is  generated  and  mixed,  the 
vaporization  time  ( tvap )  together  with  droplet  convective  time  scale 
(tconv)  are  two  important  scales  influencing  the  lift-off  location.  An 
estimation  of  tconv  was  computed  as  Zi0  mm/Uzw  Here,  Uz10  repre¬ 
sents  the  velocity  different  between  measured  droplet  (Uz)  at  the 
spray  periphery  at  z  =  1 0  mm  and  UCf.  A  characteristic  chemical 
time  scale  ( tch )  for  the  ethanol  flame  in  hot  coflow  conditions 
was  evaluated  using  calculations  of  transient  counterflow  diffusion 
flames.  The  strain  rate  was  set  to  a  representative  value  of  100  s_1. 
The  calculation  started  from  the  profile  of  a  steady  non-reaction 
flamelet  (mixing  only)  and  the  chemical  time  tch  was  defined  as 


Table  4 

Estimates  of  the  characteristics  time  scales  for  different  droplet  size  classes. 


Size-class  (urn)  Cases 

H, 

H„ 

H,„ 

j_|Aux 

ljAux 

HII 

Lop  [ms] 

10-20 

0.21 

0.31 

0.35 

0.21 

0.30 

20-30 

0.64 

0.81 

0.77 

0.63 

0.77 

30-40 

1.21 

1.48 

1.62 

1.17 

1.41 

40-50 

1.91 

2.29 

2.50 

1.85 

2.18 

tconv  [ms] 

10-20 

3.47 

3.50 

4.26 

2.65 

2.61 

20-30 

1.84 

1.80 

2.02 

1.36 

1.18 

30-40 

1.20 

1.25 

1.31 

0.96 

0.90 

40-50 

0.99 

1.02 

1.05 

0.82 

0.76 

tch  [ms] 

all 

0.31 

2.32 

4.93 

0.31 

2.32 

tvap /tconv  [~ ] 

10-20 

0.06 

0.09 

0.08 

0.08 

0.11 

20-30 

0.35 

0.45 

0.45 

0.46 

0.65 

30-40 

1.01 

1.18 

1.24 

1.22 

1.57 

40-50 

1.92 

2.25 

2.39 

2.24 

2.88 

t chi tvap  [~  1 

10-20 

1.48 

7.48 

14.09 

1.48 

7.73 

20-30 

0.48 

2.86 

5.48 

0.49 

3.01 

30-40 

0.26 

1.57 

3.04 

0.26 

1.65 

40-50 

0.16 

1.01 

1.97 

0.17 

1.06 

the  time  needed  to  reach  the  moment  of  maximal  rate  of  temper¬ 
ature  increase.  The  ChemlD  code  was  used  (CHEM1D)  using  the 
Marinov  chemical  mechanism  (Marinov,  1999).  For  the  three  cases, 
fuel  vapor  at  the  boiling  temperature  was  used  as  fuel  and  the  hot- 
diluted  coflow  as  oxidizer.  The  evolution  of  the  temperature  as 
function  of  time  at  the  stoichiometric  mixture  fraction  was  used 
to  estimate  a  characteristic  chemical  time  scale  for  each  of  the 
three  cases. 

A  comparison  between  cases  Hi,  Hu  and  Hm  shows  no  substan¬ 
tial  differences  on  the  droplets  convective  time  scale.  However,  sig¬ 
nificant  differences  are  observed  in  the  droplets  vaporization  and 
in  the  chemical  time  scales.  For  droplets  smaller  than  30  pm,  tvap 
is  lower  than  tconv  by  almost  one  order  of  magnitude.  This  indicates 
that  substantial  vaporization  of  droplets  smaller  than  30  pm 
occurs  prior  to  the  lift-off  location.  As  expected,  larger  droplets  sur¬ 
vive  the  near-injector  region  due  to  their  larger  Stokes  number  and 
spread  into  the  hot-diluted  coflow. 

For  the  highest  coflow  temperature  (case  H,),  droplets  smaller 
than  10  pm  have  tvap  of  the  same  order  of  magnitude  as  that  of 
tch.  Furthermore,  tvap  and  tch  are  considerably  smaller  than  tconv. 
This  indicates  that  chemical  reaction  engages  almost  immediately 
after  the  gaseous  fuel  is  generated  from  the  small  droplets  in  the 
near-injector  region  and  the  vapor  fuel  availability  is  the  limiting 
factor  to  form  an  ignitable  mixture.  The  lower  coflow  bulk  temper¬ 
atures  of  cases  Hn  and  Hm  yield  an  increase  of  the  tvap  and,  conse¬ 
quently,  less  fuel  vapor  is  available  prior  to  the  lift-off  location.  A 
comparison  between  the  different  droplet  size  classes  of  case  ^ 
and  Hm  show  a  reduction  of  about  60%  in  the  vaporization  time 
scales  for  droplets  in  the  size  class  of  20-30  pm  and  of  30%  for 
droplets  in  the  size  class  of  40-50  pm  (see  bold  values  in  Table  4). 
The  chemical  time  scale  increases  by  more  than  one  order  of  mag¬ 
nitude  and  is  noticeably  bigger  than  tvap  or  tconv.  It  thus  can  be 
understood  that  for  case  Hu  and  Hm,  both  tvap  and  tch  contribute 
to  an  increase  of  the  lift-off  height,  however,  tch  has  a  more  prom¬ 
inent  effect. 

The  impact  of  the  liquid  mass  flows  on  the  lift-off  height  can  be 
assessed  by  comparing  case  H,  with  Hfux  and  Hn  with  H^ux.  Increase 
of  the  liquid  mass  flows  (and,  therefore,  higher  liquid  jet  momen¬ 
tum)  yields  a  decrease  of  tconv  due  to  the  higher  droplet  initial 
velocity.  However,  only  a  slight  change  is  observed  in  the  vaporiza¬ 
tion  rates.  Hence,  the  concentration  of  fuel  vapor  in  the  droplets 
wake  until  z  =  10  mm  is  smaller  and,  consequently,  a  larger  droplet 
path  is  required  to  form  a  flammable  mixture.  This  manifests  itself 
as  an  increase  of  the  lift-off  height. 


5.  Conclusions 

A  comprehensive  study  was  carried  out  of  ethanol  sprays  issu¬ 
ing  into  hot-diluted  coflows  with  different  conditions.  LDA,  PDA 
and  CARS  techniques  were  employed  to  study  spray  flame  struc¬ 
ture  and  provide  an  explanation  of  the  stabilization  characteristics 
of  sprays  combustion  under  Mild  conditions.  An  extensive  data¬ 
base  was  obtained  for  validation  of  numerical  simulations. 

For  pressure-swirl  sprays  with  identical  Weber  number  the 
results  indicate  that  within  the  coflow  conditions  studied  the  ini¬ 
tial  droplet  velocity  and  d32  distribution  do  not  differ  significantly. 
Small  droplets  are  present  at  the  center  and  quickly  evaporate  as 
the  spray  evolves  downstream.  Large  droplets  with  relatively  high 
Stokes  number  move  outwards  and  penetrate  into  the  hot-diluted 
coflow.  The  fuel  vapor  originated  from  the  large  droplets  mixes 
with  the  hot-diluted  coflow  and  a  lean  flame-front  develops  in  a 
region  of  low  droplet  concentration  at  the  spray  outer  parts.  The 
temperature  histograms  in  this  region  show  a  considerable  amount 
of  samples  over  2000  K.  Lower  coflow  bulk  temperature  yields  a 
reduction  of  the  peak  temperatures  in  the  lean  flame-front  at  the 
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spray  outer  parts  suggesting  that  the  oxidizer  stream  temperature 
is  the  driving  parameter  with  respect  to  the  peak  temperatures. 
Due  to  the  initial  droplet  concentration  distribution,  the  gas  com¬ 
position  becomes  increasingly  rich  towards  the  spray  axis.  A  sec¬ 
ondary  flame-front  significantly  weakened  due  to  the  local 
richness  of  the  gas  composition  is  also  observed  close  to  the  spray 
center  region  for  all  the  cases  studied.  It  is  shown  that  the  observed 
increase  of  the  lift-off  height  for  lower  coflow  bulk  temperatures  is 
related  to  the  increase  of  the  droplet  vaporization  time  scale  and 
an  increase  of  the  chemical  time  scale. 
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